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Studies of Natural Electric and Magnetic Fields ' 

G. D. Garland and T. F. Webster 

(February 26, I960) 

Simultaneous measurements of short-period natural electric field variations across 
western Canada are reported. From these it is indicated thai the effect of the varying 
depth to the Precambrian rocks is the dominant factor. Analysis of the simultaneous 
magnetic and electric measurements gives a resistivity for the Precambrian basement in 
excess of 30X10 5 ohm meters. 



1. Introduction 

In considering the possibilities of earth-current 
communication, a knowledge of natural electric 
fields in the earth would appeal- to have significance 
for two reasons. First, because the natural effects 
form a background of noise above which an artificial 
signal must be detected, and secondly, because 
measurements of natural fields offer a method of 
deducing the electrical properties, in place, of por- 
tions of the earth's crust. 

It has been recognized for some time that natural 
electric currents in the earth, as detected by measure- 
ment of the tangential electric field at the surface, 
bear a relationship to magnetic disturbances. 
Recently, Hessler and Wescott [1] 2 have shown the 
variation in magnitude of the electric field as a func- 
tion of the magnetic 2£-index at College, Alaska. 
However, the pricise relationship between the surface 
electric and magnetic disturbance fields at some 
point on the earth depends on both the distribution 
of sources and on the variation of conductivity 
within the earth. It is only if these effects can be 
separated that useful information on conductivities 
can be obtained. 

An important simplification is possible if the 
source of disturbance is extended laterally, so that 
the arrival at the earth's surface may be approxi- 
mated by a plane wave. It is well known that the 
tangential electric and magnetic fields, in perpendic- 
ular azimuths, say E x and Hy, bear a simple relation- 
ship at the boundary of a uniform conductor. Then- 



EJH y =(2cT)-* 



(1) 



where a- is the conductivity and T is the period of 
the disturbance. Cagniard [2] has suggested that 
this approximation is valid for the earth, and that 
measurements of the fields will lead to a useful 
effective or apparent conductivity. However, Wait 
[3] has pointed out that departures from plane- 



1 Contribution from Department of Physics, University of Alberta, Edmonton, 
Alberta, Canada; paper presented at Conference on the Propagation of ELF 
Radio Waves, Boulder, Colo., Jan. 25, I960. The experimental work was sup- 
ported by the National Research Council of Canada. 

2 Figures in brackets indicate the literature references al the end of this paper. 



wave conditions can be treated by a series expansion 
for E x jH y . Retaining the first two terms gives 

EJHy^+rjfty'Hyimiyldy'-^Hyldx 2 

+ 2d 2 H x ld x dy] (2) 

where 7=io7xco— e/xco 2 is the intrinsic propagation con- 
stant of the earthj and rj^ifico/y is its characteristic 
impedance. Additional terms in the series involve 
higher inverse powers of 7, where I/7 is characteristi- 
cally some tens of kilometers. The simplified ex- 
pression of eq (1) will be valid if the second deriva- 
tives involved in eq (2) are sufficiently small. 

2. Regional Telluric Measurements 

It is obviously important to know the extent to 
which the fields may vary with position on the earl h 
as a result of changing source distance. In order 
to examine this, and also the effect of variations in 
earth conductivities, we have compared earth-current 
disturbances at a number of points across western 
Canada (fig. 1). The Dominion Observatory main- 
tains an earth-current and magnetic recording sta- 
tion at Meanook, and observations on the electric 
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Figure 1. Location of telluric observation points, indicated 
by solid circles, and the Meanook magnetic observatory. 
Inset, location of well where vertical currents have been detected. 
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field variations at the other points shown could be 
compared with simultaneous values at this station. 
It is the usual practice in earth-current measure- 
ments to record the changes in potential across two 
perpendicular lines (1,000 ft is a convenient length 
as the disturbances are of the order of millivolts per 
kilometer), so that any disturbance from the normal 
potential gradient may be treated vectorially. 
Within the limit of time resolution possible in these 
measurements, which were made with a paper speed 
of 4 in./hr, it was found that events could be readily 
correlated from all field points to Meanook. The 
average magnitude of a disturbance, as compared to 
that recorded at Meanook, varied from station to 
station, as shown on figure 2. These average mag- 
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Figure 2. Ratio of telluric disturbance vector at Meanook 
to that at other stations, and depth of Precambrian basement. 
Vertical lines indicate standard deviations. 



nitudes were obtained by forming the resultant of 
the two component measurements for a number of 
disturbances, of period 1 min and greater, observed 
during several hours of recording at each point. On 
the same diagram is plotted the depth to the crystal- 
line Precambrian basement, which is formed of rock 
of much higher resistivity than the overlying sedi- 
mentary rocks. There is a fairly good correlation 
between the two profiles, suggesting that the mag- 
nitude of disturbance of the electric field at any 
point is inversely proportional to basement depth. 
This would indeed be the case if source effects were 
completely uniform at all points, and if the basement 
were completely nonconducting. For then, if h is 
the basement depth, and a the uniform conductivity 
of the sedimentary rocks 



<rhE=I 



(3) 



where / is the constant total current flowing in a 
strip of unit width. At the relatively low frequencies 
considered here the skin depth is much greater than 
any basement depth encountered. It was concluded 
from these results that the major reason for the 
electric field variation with location was change 
in effective conductivity, and that source distance 
effects were minor. The approach employed here is 
an application, on a very large scale, of a method of 
geophysical exploration first described by Migaux 
[4] as the telluric method. It should be noted from 
figure 1 that the stations were located very nearly 
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Figure 3. Portion of simultaneous magnetic and telluric records. 
Note the onset of higher frequency events in hoth records at about 2.45 G.m.t. 
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along a line of constant magnetic latitude, and differ- 
ent results could well be expected for points at 
different distances from the auroral zone. 

3. Magneto-Telluric Investigations 

The apparent uniformity of the natural electric 
field suggested that Cagniard/s approach to deter- 
mine conductivity at various depths could be useful. 
Simultaneous measurements, as in figure 3, of electric 
and magnetic field variations were made at points 
where the depth to the resistant Precambrian rocks 
was known. At the time these measurements were 
made, a flux-gate magnetometer with a limit of 
sensitivity of 1 gamma (10~ 5 oersted) was used, and 
all recordings were made with pen recorders. The 
principle involved in Cagniard/s analysis is the 
extraction of a value of apparent conductivity, by 
eq (1), for a series of disturbances of different periods. 
Because of the variation in depth of penetration 
with period, the apparent conductivity determined 
will he a function of period, and the form of this 
functional relationship will depend on the variation 
of conductivity with depth. Theoretical curves of 
apparent conductivity against period have been 
computed lor cases of simple layering by Cagniard. 
It is apparent that the greatest labor involved is in 
the analysis of the records into discrete frequency 
components. 

In future observations, it is intended that this 
analysis will be made by digital recording and ma- 
chine computation. However, the values of apparent 
resistivity plotted logarithmically against period in 
figure 4 were obtained by Fourier analysis of the 
visual records. On the diagram are plotted also two 
theoretical curves for a two-layer structure. The 
upper layer is taken to be the sedimentary section, 
whose thickness is about 1.5 km at this location, with 
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Figure 4. Apparent resistivity by magneto-telluric analysis 
as a function of period. 

The curves (a) and (b) are the theoretical two-layer variations for basement 
resistivil les of infinity and 30X10 4 ohmmeters, respectively. 



a resistivity, estimated from well logs, of 30 ohm-m. 
Since the scales are logarithmic, small errors in these 
assumptions affect the position, but not the form, of 
the curves. The two theoretical curves are plotted 
for basement resistivities of 30X10 4 ohm-m and 
infinity, respectively. While there is considerable 
scatter in the observed points, the general increase of 
apparent resistivity with period is indicated. In 
particular, there is a good indication that the base- 
ment resistivity is at least 30X10 5 ohm-m. The 
observed points indicate a rate of increase of ap- 
parent resistivity greater than the theoretical curves, 
which suggests that the two-layer ease is an over- 
simplification. However, it does appear that an 
analysis of magnetic and electric fields of this type, 
especially with more adequate separation of the fre- 
quency components, will provide useful information 
on crystal conductivities. Niblett [5] has extended 
the analysis to longer periods, using records from the 
\ lea nook observatory, and has detected an increase 
in conductivity at the base of the crust. 

4. Vertical Earth Currents 

There is some interest in observing the vertical 
component of electric held within the earth, because 
this will indicate tin 4 amount of deflection of the 
horizontal currents by local variations in conduc- 
tivity, and also because of possible correlations with 
atmospheric measurements. Electrodes at depths 
of 800 and 3,800 ft have been placed in the well 
shown in figure 1. and a record of the potential 
gradient changes is available for a period of about 
1 yr. Disturbances in the vertical component 
amounting to a few percent of the horizontal electric 
field have been measured, but the ratio of vertical 
to horizontal field shows a strong dependence on the 
azimuth of the horizontal current flow. At this 
particular location, the vertical component is greater 
when the disturbance current flows north-south, as 
indicated by the plot of figure 5. This strongly 
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Ficcre 5. Hallo, obtained from well observations, of vertical 
to horizontal electric field, plotted, radially, as a function of 
azimuth of horizontal flow, for a series of pulsations from 
///,- northerly half-circle. 
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suggests lensing of the horizontal currents by a 
body of anomalous conductivity north or south of 
the observation well. However, a portion of the 
observed vertical electric field component in the 
earth may be associated with the earth-atmosphere 
current. As the normal component of current 
density is continuous across the surface, values 
of 6X10 13 ohm-m and 30 ohm-m, respectively, 
for air and earth resistivities indicate a gradient in 
the earth only 5X10 -13 times the atmospheric 
gradient. During a period of atmospheric distur- 
bance, however, the potential gradient in the air 
may be many times the fair-weather value of 100 
v/m, and the corresponding effect in the earth may 
be detectable. It is hoped that simultaneous meas- 
urements in three components, in the earth, and 
measurements in the atmosphere at this site will 
indicate the manner in which the air-earth current 
flows into, or out of, the earth. 
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